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Abstract 

We address the question of the observed pinning of i (110] ordinary 
screw dislocations in 7-TiAl which leads to the characteristic trailing 
of dipoles in the microstructure. While it has been proposed that these 
may be variously intrinsic or extrinsic in nature, we are able to rule out 
the former mechanism. We do this by means of very large scale, three 
dimensional atomistic simulations using the quantum mechanical bond 
order potential. We find that the kink-pair formation energy is large — 
6 eV, while the single kink migration energy is conversely very small — 
0.13 eV. Using these, and other atomistically derived data, we make 
kinetic Monte Carlo simulations at realistic time and length scales to 
simulate dislocation mobility as a function of stress and temperature. 
In the temperature range of the stress anomaly in 7-TiAl, we deter- 
mine whether one or several of the pinning and unzipping processes 
associated with generation of jogs are observed during our simulations. 
We conclude that the pinning of ordinary dislocations and anomalous 
mechanical behaviour in 7-TiAl must be attributed to a combination 
of extrinsic obstacles and extensive cross-slip in a crystal containing 
impurities. 

1 Introduction 

It is well known that the deformation behaviour of 7-TiAl is complex due to 
the tetragonality of the Llo structure. The 7 phase of TiAl can deform by 
several deformation modes: glide of ordinary dislocations, glide of superlat- 
tice dislocations and twinning [lj. The b = ^ (110] ordinary dislocations in 
deformed TiAl alloys exhibit a unique morphology consisting of numerous 
pinning points along the dislocation line aligned roughly along the screw dis- 
location direction, and bowed-out segments between the pinning points pQ. 
Post mortem TEM analyses [2 El H] have shown that the screw segments 
separated by pinning points are lying in different parallel (111) planes as 



a result of a cross-slip mechanism. Extensive cross-slip was also observed 
in in situ straining experiments [5]. Yield stress anomaly (increase of the 
yield stress with temperature) has been observed in 7- TiAl compounds and 
ascribed to the | (110) {111} slip system [SEIE]- There is compelling mi- 
croscopic evidence that pinning plays a prominent role in the anomalous 
mechanical behaviour. The determination of the subsequent cusp unpin- 
ning mechanism is of fundamental importance in the interpretation of the 
stress anomaly. 

Deep cusps and the trailing of dipoles suggest that the dislocations have 
many jogs. At the microscopic scale, the pinning points have been analysed 
as jogs aligned along the screw direction, whose density increases with tem- 
perature in the domain of the stress anomaly [8, 9| • The exact nature of the 
obstacles developed on these cross-kinks seems to be controversial. Viguier 
et al. [3J and Sriram et al. [I] have concluded that the pinning points are 
formed as a result of intrinsic processes involving single or double cross-slip 
mechanisms leading to formation of jogs. The three dimensional structure 
formed during collision of kinks during their lateral propagation cannot move 
in the direction of the dislocation glide. Explanations based on the activa- 
tion of cross-slip suffer however from not explaining satisfactorily a number 
of experimental results: 

1. Non-screw ordinary dislocations, Shockley dislocations involved in twin- 
ning and superlattice dislocations are also anchored at many points |1Q[ 

mi nana. 

2. Ordinary dislocations moving in the same plane are able to anchor at 
exactly the same location in the sample [5]. 

3. An annealing treatment which precipitates interstitial atoms in excess 
decreases the density of cusps and debris [2] . 

To explain these experimental results, it was concluded that pinning 
points are extrinsic in nature and due to some chemical heterogeneity [21 El 
[T5] . Messerschmidt et al. [TO] and Morris [15] have ascribed these pinning 
points to small extrinsic obstacles like oxygen atoms or AI2O3 precipitates. 

Understanding the glide mechanism of ordinary dislocations in different 
stress and temperature ranges can help us to determine which of its main 
features is likely to account for the strength anomaly of 7-TiAl. In addition, 
the behaviour of these ordinary dislocations deserves particular attention 
since they are responsible for most of the deformation in two-phase lamellar 
TiAl alloys QUI El El US] ■ 

The present work aims at exploring the mobility of ^ (110] ordinary 
screw dislocations in 7-TiAl through simulation of the specific mechanisms 
of motion of an individual dislocation in the L1q structure. 
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Due to the high Peierls barrier |16l I20j . the motion of ordinary screw 
dislocations in TiAl is believed to be controlled by nucleation, migration 
and annihilation of kinks. When the stress on a dislocation is lower than 
its Peierls stress, the dislocation line stays at rest in a given lattice position, 
interrupted by the thermally assisted process of kink pair nucleation. The 
separation of kinks under the influence of stress and thermal activation re- 
sults in translation of the entire dislocation line to the next lattice position. 
The overall dislocation movement is a cumulative effect of a large number 
of individual kink events. 

Dislocation dynamics (DD) simulations |2H [22] employ simple rules that 
dictate the motions of dislocations which are represented as interconnected 
small straight segments in an elastic continuum connected through nodes. 
Local stresses resulting from applied loading and internal stresses are com- 
puted on each of those segments. The dynamics of dislocation motion is 
thus reduced to the dynamics of the nodes prescribed by the mobility rules. 
Such a description replaces the true kink dynamics of dislocation glide with 
a mean field measure of the average mobility of a dislocation line. The re- 
liability with which the DD simulations mimic the plastic behaviour of the 
crystal depends on the accuracy of the mobility rules. 

In order to overcome these limitations and to obtain a more realistic 
description of the dynamics of a dislocation line over long time scales we 
carry out kinetic Monte Carlo (kMC) simulations of the mobility of | (110] 
dislocations. Several researchers have performed kMC studies of disloca- 
tion behaviour to provide a more realistic link between kink dynamics and 
the averaged dynamics of a dislocation line [231 Ell ES] • The advantage of 
this approach is that it replaces arbitrary assumptions about the nature of 
dislocation mobility with input based upon microscopic understanding. In 
addition, kMC can achieve extended time scales for simulation of the dis- 
location mobility. kMC simulations use atomistic results for the nature of 
the core energetics, short ranged dislocation segment interactions, and rate 
theory to simulate dislocation mobility as a function of stress and tempera- 
ture. 

In our simulations we assume that the ordinary dislocation moves by 
the kink mechanism through adjacent Peierls valleys. We examine ordi- 
nary dislocation mobility in TiAl at fixed stresses and temperatures based 
upon atomistic input. The rates of kink pair nucleation and kink migra- 
tion events are expressed in terms of their respective energy barriers, which 
are computed using atomistic simulation techniques. The kMC model al- 
lows dislocations to cross-slip onto secondary glide planes. Thus, motion 
of screw dislocations becomes three dimensional and can involve a number 
of unusual mechanisms affecting dislocation mobility. In the temperature 
range of the stress anomaly, we check whether one or several of the above 
mentioned pinning and unzipping processes associated with generation of 
jogs are observed during simulations. 
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2 kMC model 



Ordinary screw dislocations in TiAl do not dissociate into a planar structure 
due to the high energy of the complex stacking fault (CSF) in 7-TiAl |20| . 
Instead, the non-planar core of these dislocations spreads symmetrically on 
two cross-slip {111} planes (Fig. 1) and their glide is not confined to a 
single glide plane. Cross-slip should readily occur by which screw disloca- 
tions change their glide planes. In fact, each screw dislocation segment can 
nucleate kink-pairs on either of the two {111} planes. 

The kMC model does not consider any details of the core structure. It 
focuses on dislocation motion on length and time scales far greater than 
those of atomistic simulations. The key idea of kMC is to treat dislocation 
motion as a stochastic sequence of discrete rare events whose mechanisms 
and rates are computed within the framework of transition state theory. 
Dislocations are represented as interconnected small straight segments in an 
elastic continuum. Local stresses resulting from applied loading and internal 
stresses are computed on each of those segments. 

In the present model, we study the ordinary screw dislocation repre- 
sented by a piece- wise straight line stretched along the [110] direction. While 
the dislocation has, on average, a screw orientation, it consists of screw (S) 
and edge (E) dislocation segments such that kinks on the screw dislocation 
are perfect edge segments. E-segments all have the same length h, the unit 
kink height, while S-segments can be of any length. In the present work we 
choose for the length of the smallest S-segment a = b, where b = \ (110] 
is the Burgers vector of the screw dislocation. The time and 3D space is 
discretised in the form of a square grid. Here, a and h are the grid spacings 
in the screw and edge directions, respectively. 

The screw dislocation in TiAl was allowed to move on the two {111} glide 
planes that comprise the | (110) {111} slip system. Kink pairs are allowed 
to nucleate on any part of S-segments and in any of the two {111} glide 
planes intersecting the [110] direction. Once nucleated, a kink (E-segment) 
can move in its glide plane along the dislocation line until it recombines 
with another kink with the opposite sign. Periodic boundary conditions are 
applied so that kinks leaving at one end re-enter the dislocation from the 
other end. During kMC simulations, the dislocation moves under the action 
of external stress through kink-pair nucleation, migration and recombina- 
tion. This model does not allow for climb of the edge segments. 

The kMC simulation samples different classes of possible kink nucleation 
and migration events. Kink-kink annihilation is considered as a special case 
of kink migration. The kMC methods use information on the energetic 
barriers to atomic hopping or kink migration that control the motion of a 
dislocation in order to obtain the relationships between the velocity of the 
dislocation and the driving forces acting upon it. The kinetics of dislocation 
motion is completely specified by the matrix of transition rates. Transition- 
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state theory expresses the rates of kink-pair nucleation and kink migration 
events in terms of their respective energy barriers, which can be computed 
using atomistic simulations. The following form of transition-rate matrix 
element provides a physically consistent description of the energetics and 
thermally activated formation of kink-pairs [27] : 



-^k-p ( n ) = to exp 



AE n -TS- a cS bAA(n)/2 



AE n is the formation energy of a kink-pair, k-Q is the Boltzmann constant 
and co is the pre-exponential frequency factor, which is set equal to the Debye 
frequency. a e g is the local resolved stress that includes the applied stress 
and the self-stress field. The self stress is the stress exerted on a dislocation 
segment by all other dislocation segments and can be calculated from the 
gradient of the total energy of the dislocation configuration due to virtual 
displacements of the segment [27] . The total energy is computed by using 
non-singular continuum theory of dislocations [28J. The major advantage of 
this theory is that it contains no singularities. The singularity intrinsic to 
the classical continuum theory is removed by spreading the Burgers vector 
isotropically about every point on the dislocation line using a spreading 
function characterized by a single parameter r c , the spreading radius. The 
non-singular expression for dislocation energy depends on the choice of the 
cut-off parameter r c . We determined it by comparing the energy of a kink- 
pair predicted by non-singular continuum theory with atomistic data. S is 
usually taken to be 3ks |27| . The last term is the work done by the external 
applied stress and AA(n) is the change of the area of the glide plane as a 
result of kink-pair formation. 

A kink, once formed, can migrate along the dislocation line. The rate of 
migration of the kinks depends on the magnitude of the activation energy 
for kink motion W n (secondary Peierls barrier) . If this is very large the kink 
migration is thermally activated with transition rate given by 



J m = co exp 



W n - TS - a eS b 2 h/2 



(1) 



k B T 

If the secondary Peierls barrier for the kink motion is very small the kink 
migration is not thermally activated but is controlled by phonon drag. In 
such the kink velocity is proportional to the driving force expe- 

rienced by the kink. In the kMC model the thermally activated migration 
events specified by the transition rates ([T]) are replaced by continuous move- 
ment of the kinks with velocity proportional to the stress of the dislocation 
line calculated at this moment 
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2.1 Sustainable kink-pair nucleation 

The dislocations in the kMC model are discretised as sets of screw and edge 
segments. Our atomistic simulations, described in the following section, pre- 
dict that a screw segment of the dislocation line moves to an adjacent Peierls 
valley upon nucleation of kink-pair only if the distance between the kinks is 
greater than the kink width w c = 166. Only in this case can the kink-pair be 
represented by set of screw and edge segments with length h. If a kink-pair 
is formed with a width smaller than a critical width w c the height of the 
kinks is smaller than the grid spacing in the edge direction. In this case the 
kink-pair can not be represented by a continuous piecewise straight line on 
the discretised space. In addition, the elastic energy of such configurations 
can not be calculated accurately by using non-singular continuum theory of 
dislocations. A solution to this problem lies in nucleating a kink-pair only 
of finite width Wi n which is greater than the critical width w c . Such a kink- 
pair, once formed, can be represented by a continuous piecewise straight line 
on the discretised space and its elastic energy can be calculated correctly. 
To be able to do that, it is necessary to compute the rate at which such 
kink pairs should be nucleated. We determine the rate for the formation of 
a kink-pair of width w > w c by using survival probabilities analysis on the 
ID random walk. Details of this method can be found elsewhere [29J . 

The survival probability p{l —> n), which is defined as the probability 
that a kink-pair of width w = 16 expands to w = nb before it shrinks to 
w = 0, can be obtained from 

n 

p(l — s> n) = p(k — s> k + 1) 

k=l 

and p(k — > k + 1) can be calculated from the following recursive formulae 

F(k) 



p(k -> k + 1) 



1 - [1 - F(k)] p(k - 1 -t k) 



F(k) = J + {k)/ J~ (k) is the probability that upon escape from state kb the 
kink-pair will expand to width (k + 1)6. J + (k) and J~{k) are the rates of 
of forward and backward jumps from state k as follows: 



J*(k) = u exp (_W + l)-m]-TS-.*b*AWft\ {l + H(k)] 
V k BT ) 

(2) 

J-ik) = u exp ^m-^-mi-TS-^bAAjk)^ {1 + ^ 1}] 

(3) 

where H(k) = 1 if k > 1 and H(k) = if k < 0. With this definition 
the rate of embryonic kink-pair nucleation is simply J(l) = J + (0). E{k) is 



6 



the energy of two kinks separated by a distance kb. The approximation we 
employ here by describing the rates of forward and backward jumps in the 
forms d2|) and ([3|) is that the local stress remains constant over the time it 
takes for an embryonic kink-pair to expand to with nb. This approximation 
is fairly good because, under the conditions considered here, the expansion 
time is usually very small compared to the time over which the dislocation 
configuration and the local stresses change significantly. When the distance 
between the kinks kb is smaller than the critical width w c , E{k) is specified 
by the kink-pair formation energy which is calculated atomistically (sec- 
tion [321 Fig. 4, below). The rate at which kink-pairs of width Wi n b greater 
than the critical width w c appear can be obtained from the rate of embryonic 
kink-pair nucleation J(l) and the corresponding survival probability 



In our kMC model only kink-pairs with width Wi n > w c are allowed to 
nucleate on S-segments with appropriate length Lg > w c . The nucleation 
rate of such kink-pairs is given by ([!]). 

3 Atomistic simulations 

A crucial input for kMC simulations is the information from atomistic sim- 
ulations of the properties of dislocation cores and short range dislocation 
segment-segment interactions, such as the energies of a single isolated and 
kink-pair and the associated Peierls stress for kink motion. As both the 
qualitative and quantitative outcome of the kMC simulations is extremely 
sensitive to these input quantities, the accuracy of the atomistic simulations 
is of utmost importance. 

The atomistic simulation of individual dislocations requires special treat- 
ment due to the long ranged elastic field associated with them. For long 
straight dislocations, we may use a periodicity along the dislocation line, 
but should avoid using periodic boundary conditions in the other two direc- 
tions since the dislocation strain field is very long ranged. The cells used in 
these simulations, as in most real space calculations, are divided into two 
regions of relaxing and non relaxing atoms, the latter bounded by a free 
surface sufficiently far from the core to avoid image stresses. 

Pettifor's bond order potential (BOP) [301 EE] represents a numerically 
efficient scheme that works within the orthogonal tight binding approxima- 
tion with environment dependent matrix elements. The multi-atom char- 
acter of the forces, is thereby captured in a physical way that goes beyond 
the standard pair functionals (Finnis-Sinclair or embedded atom method). 
Apart from their genuine quantum-mechanical origin, BOPs have two ad- 
ditional important advantages. First, the evaluation of the energy scales 
linearly in computational time with the number of atoms, and, second, the 




(4) 
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real-space formalism avoids the need of imposing full periodic boundary con- 
ditions common to k-space methods. Both of these features are crucial for 
studies of dislocations since such simulations often require a large number of 
atoms and complex geometries. The atomistic simulations of the the prop- 
erties of ordinary screw dislocations in this study have been made using 
the BOP for 7-T1AI, constructed and extensively tested against accurate 
first-principles methods [32] . 

Molecular statics (MS) techniques, which seek to minimise the total en- 
ergy of the system at zero-temperature, can be used for relaxation of the 
inner atomistic regions of the simulation cell, where the initial atomic po- 
sitions in the computational cell are established by the conditions of linear 
anisotropic elasticity. Traditionally, fixed boundary conditions have been 
most often used in such dislocation simulations. This requires very large 
simulation cells in practice, but this method is always problematic with re- 
spect to force build-up between fixed and relaxed atomic regions. The size of 
the atomistic region can be dramatically reduced by using so-called flexible 
boundary conditions. An elastic Green's function version of such conditions 
have been developed for both 2D and 3D dislocation simulations, known as 
Green's function boundary conditions (GFBC) [33|. 134] . The GF techniques 
allow dislocation configurations to be simulated atomistically with very lit- 
tle accumulation of forces at the boundaries of the simulation cells. In this 
method, a buffer layer is introduced between the fixed outer and inner re- 
laxed atomistic regions of the simulation cell, allowing one to dynamically 
update the boundary conditions of the simulation, while dramatically re- 
ducing the size of the atomistic region. 

Using the BOP method for the description of the interatomic interac- 
tions, we have implemented GFBC to calculate the properties of ordinary 
screw dislocations through MS simulations. These properties include the 
core structure and energy, the kink-pair formation energy, the primary and 
secondary Peierls stresses as well as the kink motion activation energies. In 
this way the simulated 3D dislocation configurations can be contained in 
significantly reduced atomistic region without compromising the fidelity of 
the final core structure and energy. 

3.1 Isolated kink 

Two 2D core structures of an infinite 5 [110] screw dislocation were used to 
obtain the initial approximation to the structure of a kink. These relaxed 
configurations of the 2D ordinary screw dislocation were constructed by 
using two different elastic centres in an initial approximation to the disloca- 
tion displacement field. Anisotropic elasticity theory was used to introduce 
a 2 [110] screw dislocation in the simulation cell. A cylindrical crystal, one 
periodic unit in length along the dislocation line direction, [110], was con- 
structed in order to obtain the two symmetry-related core structures of the 
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ordinary screw dislocation. The x and y axes in the simulation cell were 
[112] and [111] respectively. The centres of the two core structures were 
introduced at two adjacent Peierls valleys with a relative displacement of h, 
the periodic length along the [112] direction. Atomistic relaxation was used 
to optimise the atomic positions in the active region. The 2D GF procedure 
was applied to relieve the boundary forces. One of the relaxed configura- 
tions obtained for the core structure of straight \ [110] screw dislocations is 
shown in a differential displacement plot in Fig. 1. 

The two 2D cores are then used to obtain the initial approximation to 
the structure of an isolated kink on an ordinary screw dislocation. A 3D 
cylindrical cell is constructed with a total length of 456, where 6 is the 
Burgers vector of the dislocation. The two halves of the simulation cell were 
made up of 22 periodic units of the two relaxed configurations of the 2D 
\ [110] screw dislocation respectively. The central cell of unit length 6 was 
defined as the average of the positions of identical atoms in the two 2D core 
structures [Sj . This overall length 456 is chosen to allow a smooth transition 
from one core structure to the next. Using a simple continuum estimate of 
the width of the kink [34(. I3"5] we obtain a kink width of 136. Therefore 456 
is a reasonable choice for the overall length of the simulation cell. 

Relaxation is carried out in the inner atomistic region and the GF tech- 
nique is used to relax the forces that develop in the buffer GF layer. The core 
structure of the central cell of the kinked dislocation line was characterised 
using differential displacement plots (Fig. 2). The difference between the 
energies of the kinked and straight dislocation gives the single-kink forma- 
tion energy. The formation energy of an isolated kink in 7-TiAl calculated 
by our atomistic simulations is 2.88 eV. 

In order to determine the width of the kink we plot the differences 
Az = [zs + 24 — (z\ + z-z)\/2 between the average ^-coordinates of the atoms 
designated numbers 1, 2 and 3, 4 in Fig. 2. This plot allows us to illustrate 
the position of the dislocation lines corresponding to the initial approxima- 
tion of the isolated kink and its relaxed core structure (Fig. 3). The width 
of the single-kink in TiAl predicted by our atomistic simulations is 166. This 
value is in good agreement with the continuum estimate of the width of the 
kink. 

3.2 Kink-pair 

The energy of two sharp kinks as a function of their separation L can be 
determined quite accurately when they are well separated |27| . The kink-pair 
can be regarded as being composed of straight line segments and the energy 
is calculated as the energy of the interactions between the line segments, 
including the self-interaction. The energy of two kinks can be separated 
into individual kink energies 2E\ and interaction energy E{. 
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E{L) = 2E k + Ei 

The energy of two mixed orientation, "oblique" kinks is expressed equiv- 
alently provided that kink-pair is separated by a distance L greater than the 
kink width w. When this condition does not hold E(L) can no longer be sep- 
arated into individual kink energies and interaction energy. The continuum 
estimation of the elastic energies of such configurations is difficult. 

Applying the BOP-GFBC technique we determine the core structure and 
energy of a kink-pair as a function of the distance between the kinks. The 
two relaxed configurations of the 2D screw dislocation were used to obtain 
initial approximations to the structure of two pure edge kinks separated by 
distances from 16 to 236. Core structures were relaxed to obtain the energy of 
the kink-pairs. The results are plotted in Fig. 4. The differences Az between 
the average z-coordinates of atoms 1, 2 and 3, 4 which illustrate the shape of 
the dislocation lines for several kink-pairs are shown in Fig. 5. The atomistic 
simulations predict that a screw segment of the dislocation line between both 
kinks transfers to the adjacent Peierls valley only if the distance between the 
kinks L is greater than the kink width, (i. e. L > 166). The formation energy 
of two kinks in this case can be determined quite accurately by using non- 
singular continuum theory of dislocations [27]. By comparing the energies 
of two kinks separated by distances L > 166, predicted by non-singular 
continuum theory, with the corresponding atomistic data we can determine 
the spreading radius (see section [2]) or cut-off parameter. The value of the 
cut-off parameter obtained in this way is r c = 0.1256. 

We made continuum calculations of the formation energy both for pure 
edge, or "perfect", and oblique kinks. In both cases we obtained similar 
results for the energy of the kinks. The relative discrepancy between the 
formation energies of perfect and oblique isolated kinks is 0.25. For kink- 
pairs separated by distance L > w we find that the relative discrepancies 
between the energies of perfect and oblique kinks is about 0.5. Taking 
into account these calculations, we approximate the real oblique kinks on 
ordinary screw dislocation in our kMC model by perfect edge segments. The 
approximation of perfect kinks is appropriate when the kink width is smaller 
than the distance between two kinks in a kink-pair. In this case, calculation 
within the non-singular continuum theory and perfect kinks approximation 
provides a correct result for the self-stress exerted on a dislocation segment. 

3.3 Kink migration 

In order to estimate the energy barrier for kink migration we started with 
a simulation block containing the fully relaxed core structure of an isolated 
kink on an ordinary screw dislocation. To this we applied an incrementally 
increasing pure shear stress such that the force on the dislocation pushes 
the kink along the dislocation line. In practice this stress is imposed by 
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applying the appropriate homogeneous shear strain, which is evaluated with 
anisotropic elasticity theory. This strain was superimposed on the disloca- 
tion displacement field for all the atoms. Relaxation was carried out at every 
incremental step in the applied stress, until the kink migrates one crystal 
period along the dislocation line. At each step we calculated the change 
in the total energy of the cell as the kink migrates from its initial to final 
position. In particular, the maximum energy difference along the path is the 
energy barrier of migration of a single kink. Fig. 6 shows the energy of an 
isolated kink during its motion between two adjacent Peierls valleys. The 
energy is plotted at the nine incrementally increasing shear stresses that we 
have been working with. We obtain 0.13 eV for the energy of migration of 
a single kink on a screw ordinary dislocation in 7-T1AI. 

In contrast to high Peierls barrier for kink-pair nucleation (~6 eV, see 
Fig. 4), kink motion on the ordinary screw dislocation in TiAl should be very 
easy, as manifested in the low kink (or secondary) Peierls barrier. A high 
mobility of nucleated kinks should be expected given that the applied stress 
in our simulations (rather less than the Peierls stress ap = 0.04Caa [20]) far 
exceeds the secondary Peierls stress (O.OO23C44) determined in our atom- 
istic simulations. The low energy barrier and secondary Peierls stress mean 
that kink migration along the ordinary screw dislocation is not thermally 
activated but is limited by phonon drag. Therefore, the kink velocity in the 
present kMC simulations is calculated as 

b 

where a g is the glide component of the resolved shear stress, b the Burgers 
vector and B the phonon drag coefficient. In our simulations we use the same 
numerical value of the drag coefficient for kink motion, B = 1.5 x 10 -5 Pa s, 
as in Fivel et al. |36j. 

4 Results 

In this section we present the results of kMC simulations of ordinary screw 
dislocation motion under 0.75<7p shear stress with different orientations of 
the maximal resolved shear stress (MRSS) plane. The total length of the S- 
segments are 40006 ~ 1.1 //m while each E-segment is only 6.5 A. Since the 
typical measured distance separating two pinning points is 0.1-0.15 fim |13j . 
the choice of this length of the dislocation line in our simulations is to 
maintain enough space for intrinsic pinning of the ordinary dislocation. Two 
different shear stress orientations of the MRSS plane have been tested in 
our kMC simulations: (i) single slip conditions, and (ii) non symmetrical 
double slip conditions, each of them at different temperatures between 600 
and 1200 K. The choice of such a moderately high stress and temperature 
values in these simulations is to maintain a sufficient number of kinks on 
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the dislocation. In the case of single slip conditions the MRSS plane is 
parallel to the primary glide plane (111). The resolved shear stress on the 
primary glide plane is much larger than that in the secondary plane, and 
kink-pairs mostly nucleate on this plane. If the pinning points are supposed 
to be intrinsic, they are initiated by the collision of superkinks nucleated 
in both primary and cross-slip planes, and transformed into jogs by double 
cross-slip. In the case of non symmetrical double slip conditions the MRSS 
plane bisects two glide planes (111) and (111), making kink-pair nucleation 
almost equally probable on both planes. We expect that collision of kinks 
nucleated in both primary and cross-slip planes resulting in self-pinning 
should be readily observed in this case. 

However, for both shear stress orientations of the MRSS plane the sim- 
ulation results were very surprising, in that there is always only one kink- 
pair along the entire dislocation line. Kink nucleation is rare and existing 
kinks tend to glide along the dislocation line through the periodic bound- 
ary conditions and recombine before the next kink nucleation event occurs. 
The averaged kink-pair nucleation rate per lattice site on the glide plane is 
Jk_ p = 1.5 x 10 5 b^ 1 s _1 , and the average kink velocity is ~ 10 14 6 s _1 . 
The interplay between the rate of kink-pair nucleation and the kink mi- 
gration velocity gives rise to a particular length scale L c such that when 
a dislocation is shorter than L c , there is always only one kink-pair along 
the entire dislocation line. For typical kink-pair nucleation rate and kink 
migration velocity in our simulations we find that this length scale is more 
than 10/i m. Comparison with the typical observed distance separating two 
pinning points indicates that the probability for self-pinning and an intrin- 
sic nature of the pinning points of screw ordinary dislocation in 7-T1AI is 
negligible. 

In order to demonstrate self-pinning of the ordinary dislocation we arti- 
ficially reduced the kink migration velocity by increasing the phonon drag 
coefficient. With decreasing kink migration velocity the number of kink- 
pairs nucleated on both glide planes along the dislocation line steadily in- 
creases. At B = 5 x 10~ 3 Pa s, corresponding to a kink migration velocity 

?3 3 x 10 11 6 s _1 a few pinning points are generated along the dislocation 
line, leading to cusp formation, see Fig. 7. The rate of self-pinning and cusp 
formation along the screw direction is higher in the case of non symmetrical 
double slip conditions. The process of intrinsic pinning is initiated by col- 
lision of two kink-pairs formed in both primary and cross-slip planes. Such 
elementary jogs grow in size when more kinks pile-up on either sides of the 
initial pinning point forming superjogs. The kinks in the pile-ups on two 
sides of the pinning point belong to different parallel primary glide planes. 
If the dislocation lines on two sides of the pinning point cross each other 
they may reconnect by recombination of kink-pairs. As a result, the dislo- 
cation line is now reduced in size, leaving behind a prismatic loop as shown 
in Fig. 7. The cusps are easily released by this unzipping mechanism and 
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straight dislocations along the screw direction are restored. 

The separation between pinning point and the unzipping are functions 
of temperature and shear stress orientations of the MRSS plane. In the case 
of non symmetrical double slip conditions pinning, unzipping and formation 
of prismatic loops aligned along the screw direction is frequently observed 
in our kMC simulations. At lower temperatures and single slip conditions 
the rate of self-pinning and unzipping is less. Due to the lower unzipping 
rate, in this case we observed formation of dipoles as a result of propagation 
of screw dislocation line on both sides of the pinning points. 

5 Discussion 

The results of our simulations indicate that in the unit processes of ordinary 
dislocation glide described by the kink pair formation and migration model, 
nucleation of sustainable kink-pairs is rate limiting, and hence plays an 
important role in determining both the velocity of the dislocation and its 
glide mechanism. 

The specific ordinary dislocation mechanisms expected in the LIq struc- 
ture were explored by Fivel et al. [36J through an adaptation of a 3D meso- 
scopic dislocation simulation initially developed by Kubin et al. [37]. Using 
simple rules for kink-pair nucleation and propagation on ordinary screw dis- 
locations, the simulation reproduced the main features of the "local pinning 
unzipping" (LPU) mechanism [38j (pinning on cross-slip generated jogs, cusp 
unzipping, formation of aligned prismatic loops). It was concluded that ex- 
trinsic obstacles do not seem to be necessary to account for both the observed 
microstructures nor the stress anomaly [36] . This conclusion concerning the 
dislocation glide mechanism was a result of the the high kink-pair nucle- 
ation rate employed by Fivel et al. [36]. The high probability of thermally 
activated kink-pair nucleation followed from low kink-pair activation energy, 
0.45 eV, which they used. Since no experimental measurement of the acti- 
vation energy is available owing to the stress anomaly, the activation energy 
in 36j is simply taken to be in the range of the analytically calculated activa- 
tion enthalpy for the cross-slipping back from the secondary (111) cross-slip 
plane [9]. 

On the other hand, our atomistic simulations show that the kink-pair 
activation energy is one order of magnitude higher, which in its turn leads to 
a much lower thermally activated kink-pair nucleation rate. The interplay 
between the rate of kink-pair nucleation and the high kink migration veloc- 
ity due to low secondary Peierls barrier leads to low density of kink-pairs in 
both {111} glide planes at length scales at which pinning of the dislocation 
line should be expected. In the absence of obstacles, kinks on the ordinary 
screw dislocation have very high mobilities, limited predominantly by lattice 
drag resistance. In the presence of solute atoms, the rate of kink migration 
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is determined by the solute-kink interaction. The solute-dislocation inter- 
action can be divided into two parts depending on the distance of the solute 
from the dislocation segment; one associated with the long ranged elastic 
fields of the solute and the other with the short ranged solute-core interac- 
tion. If the solute is close to the dislocation core, the core-solute interaction 
provides a barrier Ef> to the motion of the kink. Usually kinks require ther- 
mal activation to overcome the local energy barrier E^ associated with solute 
atoms near the core. In this case, the rate of kink migration is determined 
by the magnitude of the solute-kink interaction energy E^. 

Owing to high kink-pair activation energy and high kink mobility in the 
absence of interaction with extrinsic obstacles predicted by the present atom- 
istic simulations, the probability for self-pinning is negligible, and therefore 
the nature of the pinning points on screw ordinary dislocations in pure LIq 
crystal cannot be intrinsic. It follows that the observed dislocation pinning 
must be attributable to extrinsic obstacles. These obstacles presumably re- 
duce kink migration velocity, thus increasing the density of kink-pairs in 
both {111} glide planes at length scales at which pinning of the dislocation 
line should be expected. The increased density of kink-pairs formed in both 
primary and cross-slip planes leads to higher probability for collision of two 
kinks and formation of jogs. 

The results of our simulations indicate that the pinning on screw or- 
dinary dislocations can be attributed to the relation between extrinsic ob- 
stacles and extensive cross-slip. As shown in the present study, the lateral 
propagation of kinks occurs very quickly and this implies a very small prob- 
ability of movements in the adjacent slip plane of the screw segments which 
are lying in the sessile configuration before this lateral propagation. If a 
part of a dislocation moving without any interaction with extrinsic obsta- 
cles makes a double cross-slip, the whole dislocation moves onto a primary 
glide plane parallel to the initial one since no obstacles oppose the lateral 
propagation of the edge segments lying in the cross-slip plane. If the dis- 
location moves in a crystal containing extrinsic obstacles, the number of 
kink-pairs nucleated on both primary and cross-slip planes along the dislo- 
cation line steadily increases as a result of the reduced mobility of kinks due 
to the solute-dislocation interaction. The increased density of kink-pairs 
nucleated in both glide planes leads to higher probability for collision of 
two kinks and formation of jogs by double cross-slip. The pinning points 
and dipoles dragged by screw dislocations are thus probably formed through 
this combination of extrinsic obstacles and cross-slip in a crystal containing 
impurities. 
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6 Conclusions 



In our kMC simulations of ordinary dislocation mobility in TiAl we exam- 
ined an isolated dislocation moving at applied stress by the kink propagation 
mechanism through adjacent Peierls valleys in the Ll$ crystal structure of 
7-TiAl. The rates of kink nucleation and migration events were expressed in 
terms of their respective energy barriers, which were computed using fully 
three dimensional atomistic simulations employing Green's function bound- 
ary conditions and a quantum mechanical prescription for the interatomic 
forces. The results of the present simulations indicate that the interplay be- 
tween the rate of kink-pair nucleation and the kink migration velocity leads 
to a surprisingly low density of kink-pairs in both glide planes at length scales 
at which intrinsic pinning of the dislocation line should be expected. There- 
fore, the probability of intrinsic pinning, which is initiated by the collision 
of kinks nucleated in both primary and cross-slip planes, and transformed 
into jogs by double cross-slip is negligible. 

If we artificially reduce the kink migration velocity the number of kink- 
pairs nucleated on both glide planes along the dislocation line steadily in- 
creases. The increased density of kink-pairs formed in both primary and 
cross-slip planes leads to higher probability for collision of two kinks and 
formation of jogs, demonstrating that our method is capable of producing 
the intrinsic self-pinning, albeit at unphysical values of the parameters. 

We conclude that the pinning points are not intrinsic in nature. The 
pinning of ordinary dislocations and anomalous mechanical behaviour in 7- 
TiAl must therefore be correlated to extrinsic obstacles. Amongst the other 
possible mechanisms, obstacles could reduce kink migration velocity, thus 
increasing the density of kink-pairs in both glide planes at length scales at 
which pinning of the dislocation line should be expected. 
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Fig. 1. Differential displacement plot showing the screw components 
of a i [110] screw dislocation. The arrows joining atom pairs indicate the 
relative displacements of the two atoms in the direction normal to the page; 
hence the arrows show the screw component of the relative displacement. 
The length of the arrow (in arbitrary units) is proportional to the amount 
of relative displacement. 

Fig. 2. Differential displacement plot showing the core structure of the 
central cell of the kinked dislocation line. The arrows are as in Fig. 1. This 
shows the plane of atoms perpendicular to the Burgers vector in the region 
separating two long screw dislocations having elastic centres displaced by the 
kink height, h. Hence a kink is lying along the [112] direction in the plane of 
the diagram. The core of the two \ [110] screw dislocations spreading onto 
two adjacent (111) and two adjacent (111) planes can thereby be clearly 
seen. 

Fig. 3. Core spreading of an isolated kink. The relative displacement, 
Az = [zs + Z4 — (z\ + Z2)]/2 between the average z-coordinates of the atoms 
designated numbers 1, 2 and 3, 4 in Fig. 2 is plotted against distance from 
the kink centre in units of b. In the ideal, unrelaxed case this amounts to a 
"perfect" kink. The relaxed structure represents what in the text we refer 
to as an "oblique" kink, or non pure edge kink. We see that the width of a 
relaxed kink is on the order of ten Burgers vectors. 

Fig. 4. Kink-pair formation energy as a function of the separation 
between the kinks. Here we constrain the kink-pair separation in the 3D 
simulation as described in the text and show the energy of the corresponding 
relaxed kink-pair. 

Fig. 5. The shape of the dislocation lines for kink-pairs with different 
distances w between the kinks. As in Fig. 3, we show the relative displace- 
ment, Az, this time at the cores of a kink-pair, constrained as in Fig. 4 to 
a separation w/b, this being indicated on each set of data. The dotted line 
at w = 21b refers to an "ideal", unrelaxed kink-pair. This illustrates the 
point we make in the text that for a kink-pair to be identifiable in a kMC 
simulation its separation must be greater than ~ 166, so that it extend into 
the next Peierls valley. 

Fig. 6. The energy of an isolated kink during transition between two 
adjacent secondary Peierls valleys. The applied stress is in units of the shear 
elastic constant, C44. 

Fig. 7. A typical outcome of the kMC simulations, showing the cusped 
character of the dislocation line observed at reduced kink mobility. The 
MRSS plane bisects two glide planes (111) and (111) making kink-pair nu- 
cleation almost equally probable on both planes. The process of self-pinning 
is initiated by the collision of kinks nucleated in both primary and cross-slip 
planes. The cusps are released by recombination of the kink-pairs on both 
sides of a pinning point, leaving behind debris (prismatic loops). 
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